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Energy and momentum of the elementary excitations become independent vari- 
ables in medium: energy and momentum statistical distributions are not identical. 
The momentum distribution and not the energy distribution is relevant for barrier 
penetration. The deviations of the momentum distribution from the Maxwell- 
Boltzmann energy distribution can be expressed in terms of the imaginary part 
of the self-energy of the quasi-particle. It is possible to obtain an effective Tsal- 
lis' distribution for the kinetic energy. These effects are different from static or 
dynamical screening and can have important consequences for reaction rates in 
stars. 

1 Introduction 

Theoretical calculations of nuclear rates as function of temperature, density 
and composition are fundamental ingredients of our understanding of stellar 
structure. Solar standard models (SSMs) are based on rates of the nuclear re- 
actions inside the Sun calculated,according to what is the actual experimental 
and theoretical state-of-the-art El. 

The solar structure is quite robust and even large changes of nuclear 
rates yield non standard solar models whose structure is very similar to the 
standard one; only recent precise measurements of helioseismic frequencies 
can discriminate between standard and non standard models u. In addition, 
there exist changes in chemical composition, such as those produced by He 3 
mixing □, that cannot be detected even by seismic measurements. 

The flux and energy spectrum of neutrinos from the Sun are observables 
quite sensitive to some of the reaction rates, but the possibility of oscillations 
weakens the link between the observed fluxes and the rates inside the Sun 

It is, therefore, of the outmost importance that nuclear rates inside the 
Sun, and more in general inside the stars, be accurately calculated and all 
possible effects be taken into account. 

Cross sections are fundamental for the determination of the rates: the 
astrophysical S-factors used in stellar model calculations are often extrapola- 
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tions from experimental data at higher energies than those relevant for stellar 
interiors; some reactions, e.g., the p + p reaction, cannot be measured. Al- 
though the underlying theory is generally robust, new calculations are still 
necessary. n 

As an example, a recent and more accurate calculation a of the He 3 + p 
reaction predicts a S-factor five time larger than the one used in the SSM. 
Even if this reaction remains not important for the basic structure of the Sun, 
it becomes not negligible for the high-energy part of the neutrino spectrum. In 
particular, it could be important in the interpretation of the possible excess of 
solar-neutrino events above 13 MeV, even if the latest experimental data seem 
to suggest that the excess was a statistical fluctuation. Nevertheless, there 
exist a vacuum oscillation solution to solar neutrino problem that reproduces 
the apparent seasonal variation of the temporal series of the GALLEX and 
Homestake data and that would imply a high energy distortion of the spectrum 
of the recoil electrons 13. The most distinct signature of this solution is a semi- 
annual seasonal variation of the 7 Be neutrino flux with maximal amplitude, 
future detectors (BOREXIND, LENS and probably GNO) will be able to test 
it: the accurate calculation lI of the He 3 + p reaction will be very important 
also for the theoretical interpretation of these results. 

However, nuclear rates depend not only on the reaction cross sections but 
also on the properties of the nuclear plasma. Nuclear plasmas, i.e., neutral 
systems of charged particles (ions and electrons), are complicated many-body 
systems: charged are screened, but the screening depends in general on the 
energy of the particle and the remaining interaction is long range and non 
local. The resulting spatial and temporal correlations between ions have large 
effects especially on fusion reactions, which occur between those high-energy 
ions that can tunnel the Coulomb barrier. 

An accurate theoretical determination of these rates requires a good un- 
derstanding of all possible plasma effects; at the same time, reactions that 
select ions in the tail of the energy distribution are probes of the dynamics of 
the plasma itself. 

In this contribution we only want to discuss one such plasma effect: the 
possibility that the effective momentum distribution of the fusing ions could 
deviate from the standard Maxwell- Boltzmann (MB) distribution. In particu- 
lar, we want to report on our work in progress that tries and link the deviation 
from the MB distribution to the imaginary part of the self-energy of the quasi- 
particle states. 

While this effect is potentially important for many astrophysical phenom- 
ena (the solar neutrino problem, brown dwarves, dark matter distribution), 
we shall focus on parameters relevant to the solar core. 



Submitted to World Scientific on November 30, 2000 



2 Non ideal plasmas 



The solar core is a weakly-nonideal plasma where: 1) the mean Coulomb 
energy potential is of the order of the thermal kinetic energy; 2) the Debye 
screening length Rrj ss a (interparticle distance): Debye- Hiickel conditions 
are only approximately verified; 3) it is not possible to separate individual 
and collective degrees of freedom; 4) the inverse solar plasma frequency (t p i = 
oj~l — \J m/Aime 2 ~ 10~ 17 ) is of the same order of magnitude of the collision 
time t co u = / _1 = (no~v); 5) particles loose memory of the initial state only 
after many collisions: the scattering process cannot be considered Markovian; 
6) the time needed to build up again the screening, after hard collisions, is 
not negligible ml 

At the thermal equilibrium reacting ions are usually described as quasi 
free particles with Maxwell-Boltzmann (MB) velocity distribution. But many- 
body effects inside the plasma could cause deviations from a pure Maxwell- 
Boltzmann statistics for the effective degrees of freedom. Because reacting 
ions belong to the high momentum tail of the distribution, at least for fusion 
reactions between charged ions, even tiny deviation from the MB tail can 
cause large modifications (enhancement or depletion) of the rates. 

The value of the collision frequency / determines the possibility of two 
different effects that produce important deviations from the Maxwellian dis- 
tribution Fm{p) at high momenta: 

Q) Quantum uncertainty effect. 

When the Coulomb collisional frequency is large (hf > kT) the ions 
cannot be considered as quasi-free particles: the energy and momentum 
distributions are different and one must decide which one is relevant for 
the reaction rates. The fact that the two distributions are not equivalent 
is related to the finite life-time of the quasi-particles and to the quantum 
uncertainty. Since nuclear rates should be evaluated averaging the quasi- 
classical cross section o~{p) over the momentum distribution, rather than 
the energy distribution, even if the energy distribution is -Maxwellian, the 
effective distribution can acquire a non-Maxwellian tail 

q) Weak nonextensiyity effect. 

Tsallis statistics 113 with entropic parameter q can describe systems that 
are not extensive due to long-range interactions or non-Markovian mem- 
ory effects; the energy distribution itself deviate from the standard free- 
particle statistics. When deviations are small (q sa f ) the correction (en- 
hanced or depleted tail) can be described by the factor exp — {vr) 
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Deviations from the Maxwellian tail due to either Q or q effect (or both) 
may lead to strong increase or decrease of the nuclear rates in the solar core 
(non standard solar models due to large changes of the nuclear rates and their 
implications for the solar neutrino problem are described in Ref. □). 

In this contribution we discuss only the Q effect, i.e., the deviation of the 
momentum distribution relative to the energy distribution. As we shall see, 
this effect leads only to an enhanced tail. 

3 Quasi-particle momentum distribution 

Many properties of interacting systems can often be described by weakly in- 
teracting excitations or quasi-particles. The energy-momentum dispersion 
relation (position of the pole of the one-particle Green's function) of these 
excitations is found by solving 



where Y,(ui,p 2 ) = E# + zEj is the self-energy of the one-particle propagator. 

In the approximation of a constant real part of the self-energy E#, we ob- 
tain the shift of the energy due to the static mean-field (in plasma it produces 
non-dynamical screening) . 

The p 2 dependence of E# reflects the spatial nonlocality of the effective in- 
teraction and may be understood qualitatively by considering the nonlocality 
of the exchange term of the Hartree-Fock potential, while the to dependence 
reflects the nonlocality of E in time. 

As long as the imaginary part Ej is zero, there exist a one-to-one cor- 
respondence between the energy to and momentum p (or kinetic energy 
p 2 /(2m)). In real systems the imaginary part of the self-energy of the quasi- 
particle is not zero and energy and momentum become independent variables; 
however, they are still strongly correlated when Ej is small (only if Ej is small 
the concept of quasi-particle is useful) . 

Barrier penetration is determined by the momentum of the (quasi-parti- 
cle and not by its energy (when they do not coincide). 

In this preliminary presentation, we shall restrict ourselves to the case of 
an energy distribution that is Maxwellian: P(E) ~ exp(— j3E). 

If we are given the relation between E and p 2 in the form F(E,p 2 ), the 
momentum distribution is obtained 




(1) 




(2) 
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For free particles F(E,p 2 ) — S (E — p 2 /(2m)) and, therefore, 

P(p 2 ) = e- (i ^ . (3) 

If S/ = and T, R is constant, then F(E,p 2 ) = 6(E-p 2 /(2m) - T, R ) 
and the distribution is still Maxwellian 

P(p 2 ) = e \ J . (4) 

If S/j is not constant but can be expanded in the region of interest 

dp 2 iP Po) + duj 

the result is 



X r (cu,p 2 )=V r + ^-(p 2 -p 2 ) + ^(lo-lo ) (5) 



P(p 2 )=e V- S , (6) 
where 

m . = ra ( 1+2m | f )-'( 1 _^) , (T) 

In all the above cases the distribution for the variable p 2 is still Maxwellian 
(in general it follows the energy distribution), even if the shift in energy or the 
effective mass can have important phenomenological consequences (screening, 
level densities, etc.). 

When the imaginary part of the self-energy cannot be disregarded 
(£/ > 0), there appear deviations from the Maxwellian distribution. For 
the sake of discussion let us consider the following relation between ui and 
e p = p 2 /(2m*) + 

where the parameter v characterizes the deviation from a (^-function. This 
function is normalized 

/ dEF„(E,p 2 ) = —— / dxx 1 '"-^-* = 1 , (9) 
Jo 1 w^j Jo 



and in the limit v — > + 

log(^) = i 



, E i E 

1 + log — 



-Uog(2nvE 2 ) + 0(v) (10) 
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The function between bracket f(x) — 1 — x + log(:r) < and is zero only for 
x = 1; therefore, if E ^ e p , lim„^ + log(.F„) = — oo and lim„_, + F v = 0. If 
E = e p , lim„^ + F v = ^^ nvE — °°- In summary 

lim F„(E,p 2 ) =d(E-e p ) . (11) 

If v is small but not zero, we try the Ansatz E/e p = 1 + a^Jv into the equa- 
tion ([!(]) and find 

log(F„) =-\~\ \og(2nvE 2 ) + O(iz) (12) 

where the terms that have been dropped are really of order v as long as 
a = (E/e p — V)l\Jv remains of order one, i.e., as long as E/e p — 1 does not 
becomes large compared to ^fv (remains of order \fv). 

From the definition a? = (Ejt p — l) 2 /^ we can rewrite Eq. ( |l2| ) 

• (13) 



Since \fvt p is the width of the distribution and the imaginary part of the 
self-energy £/ is also proportional to the width of the distribution (when it is 
in Lorenzian form), it is plausible that 

„- C {^) . (14, 

at least in the limit £/ <C e p ; the constant C depends on the precise defini- 
tion of the limit of the quasi-particle (we are working on the microscopical 
derivation of such kind of relation and of C) . 

If we instead calculate the distribution of e p from the relation F u {uj,p 2 ) 
between energy and momentum for general v (without expansion for small v) , 
we find 




r°° r°° 11/ 

/ du;e-^F(cj,p 2 ) = / dcje~^-—- 

Jo Jo ruH w V 



1 1 / UJ 

ve p 

l + vfle,)-* , (15) 



which is the Tsallis' distribution with v = q — 1 > 0. 

It is not possible to obtain a Tsallis' distribution with v = q— l<0by this 
kind of effect, since the broadening of dispersion relation between energy and 
momentum physically has the effect of increasing the tail of the distribution: 
it cannot "cut" the tail. 
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At least in the limit of small £/, we have given a phenomenological inter- 
pretation of the parameter q in the Tsallis' distribution 

,-^c(|) J . (le, 

4 Width comparison 

We are studying a more general approach that would allow to compare dif- 
ferent relations between p 2 and E, at least in the asymptotic limit that each 
relation F(E,p 2 ) tends to a ^-function, F(E,p 2 ) — > S(E — e p ). As long as 
we consider one-parameter generalizations of the 5-function, there should be 
some mapping from one parameterization to the other. We have already seen 
how the relation in Eq. (||) and a Gaussian relation coincide in the limit of 
small width. 

An other interesting case is the Lorenzian, which also becomes a 8- 
function in the limit of vanishing width. In fact, the quasi-particles dispcpion 
relation becomes in dense media (at least in the limit of large life-time) c3 

<w) = ~— — ; — : ^ 2 i — : ^ » ( 17 ) 



7r [(e - e p - A(e,e p )) 2 + g(e. 



where e p — p 2 /2m, A(e, e p ) and g(e,p) are the real and imaginary parts of the 
one-particle retarded Green's function self-energy. 

For weakly non ideal plasmas one can show that A w kTT/2, where T is 
the plasma parameter (r = e 2 /RokT) and g oc hf. At non zero value of g, 
a nonexponential tail appears, in. the distribution function Fq (p) . For large 
momenta, it has been found liSEd that 

+ 27r s 

where /i is the chemical potential and Fm{p) the Maxwellian distribution. 



F Q {p)^F M {p)+ } g h ^e^^ (18) 



5 Conclusions 

We have shown that even when ions have a MB energy distribution the finite 
life-time of the quasi-particles in the plasma can produce a non-Maxwellian 
momentum distribution. Since the tunneling probability between charged 
ions must be evaluated using the momentum distribution, the reaction rates 
are effectively obtained by using distributions that depart from the MB one. 
This departure from the MB distributions can be calculated if one knows the 
spectral dispersion relation between energy and momentum. 
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In particular, we have shown that it is possible to have relations between 
energy and momentum that yield momentum distributions of Tsallis type 
with entropic parameter q > 1, which corresponds to an enhanced tail. 

In this framework, we have suggested a possible interpretation of the 
parameter q in terms of the imaginary part of the quasi-particle self-energy, 
at least in the limit of q — ► 1 (small deviations from MB). 
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